Comparison of a Converted Helmholtz and 
Euler Model for Impedance Eduction in Flow 

Willie R. Watson* and Michael G. Joncs t 
NASA Langley Research Center, Hampton, Virginia 23681-2199, USA 


Impedances educed from a well-tested convected Helmholtz model are compared to that of a recently de- 
veloped linearized Euler model using two ceramic test liners under the assumed conditions of uniform flow 
and a plane wave source. The convected 1 Ielmholtz model is restricted to uniform mean flow whereas the lin- 
earized Euler model can account for the effects of the shear layer. Test data to educe the impedance is acquired 
from measurements obtained in the NASA Langley Research Center Grazing Incidence Tube for mean flow 
Mach numbers ranging from 0.0 to 0.5 and source frequencies ranging from 0.5 kHz to 3.0 kl Iz. The unknown 
impedance of the liner is educed by judiciously choosing the impedance via an optimization method to match 
the measured acoustic pressure on the wall opposite the test liner. Results are presented on four spatial grids 
using three different optimization methods (contour deformation, Davidon-Flelcher Powell, and the Genetic 
Algorithm). All three optimization methods converge to the same impedance when used with the same model 
and to nearly identical impedances when used on different models. An anomaly was observed only at 0.5 kHz 
for high mean Mow speeds. This anomaly is likely due to the use of measured data in a flow regime where 
shear layer effects are important but are neglected in the math models. Consistency between the impedances 
educed using the two models provides confidence that the linearized Euler model is ready for application to 
more realistic llows, such as those containing shear layers. 
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= finite element system matrix mid vector of source effects. Pa 
= diameter and depth of the parallel cylindrical channels composing the lest liner, m 
= ambient sound speed, m/s, and free space wave number, 1 / m 
= wall error function, Pa 2 

= source frequency, Hz, and angular frequency, s -1 
= height, width, mid length of duct, m 
= y/— 1, unit imaginary number 

= location of the leading mid trailing edge of wall lining, m 
= mean flow Mach number 

= fully developed mean flow Mach number mid centerline Mach number 
= average mean flow Mach number over three axial stations 
= order of [A] mid number of lower wall measurement points 
= mean flow pressure. Pa, temperature, K, .and density, kg/m 3 
= frequency domain acoustic pressure. Pa 

= dimensionless acoustic resistance and reactance normalized by poCo 
= time, s 

= frequency domain axial mid transverse acoustic particle velocity, m/s 
= transverse, spanwise, mid axial coordinate, m 
= fl 4- i%, dimensionless acoustic impedance normalized by pof'o 
- vector of unknown node pressures, Pa 
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Subscripts: 

0 

exit 

FEM, Meas 
/ 

ref 

s 

c 


= mean flow quantity 
= exit plane quanity 

= finite element computed and measured acoustic pressure 
= index of summation 
= standard day condition at 295 K 
= source plane 
= centerline 


Superscripts: 

* = complex conjugate 


I. Introduction 

E FFICIENT duct liners for broadband acoustic noise suppression remain critical to the development of environmen- 
tally acceptable commercial aircraft. Previous experiments 1 have shown that significant changes in duct liner 
impedance can occur as a result of grazing flow. Because of the need to better understand this phenomenon, NASA 
Langley Research Center (I^iRC) has invested a considerable effort into the development of convcctcd Helmholtz 
equation (CHE) 2-4 and linearized Euler (LEE) 5 models to educe impedance in grazing flow. Although the CHE model 
is restricted to uniform flow, the LEE model can account for the effects of the shear layer. In both models, the test 
section in the InRC 51x51 mm Grazing Incidence Tube (GIT) is used to acquire the acoustic and mean flow data 
needed to educe the impedance of the test liner. Impedances educed using the CHE model have been extensively 
reported 2-4 and compare favorably to normal incidence impedance tube results at Mach zero. Inconsistencies in the 
impedance educed from the CHE model have been observed at the lowest frequency (i.e., 0.5 kHz) at high mean flow 
speeds. To date, no satisfactory explanation for this anomaly has been given. 

In 2002, results from the CHE model were compared to results of methods used in various IPS. aeroacoustics 
laboratories. 6 Impedance data comparisons acquired from this multi-laboratory study suggested the need to incorpo- 
rate the shear layer and 3-D acroacoustic effects into the impedance eduction methodologies. Subsequently in 2003, a 
95-microphone window and improved mean flow' suction devices w'ere installed in the GIT 7 to obtain acroacoustic and 
improved mean flow measurements for impedance eduction in 3-D sound fields with shear layers. In 2005, a model 
w'as developed 4 for educing liner impedance using acroacoustic and mean flow' data obtained from the modified GIT. 

Although the LEE model 5 can account for the shear layer, it has yet to be extensively validated even in uniform 
mean flow. One would expect that impedances educed using the CHE and LEE models w ould be identical for the same 
frequency, source condition, mean flow field, and liner. In addition, several time domain computational aeroacoustics 
(CAA) models that are being developed to model the effects of duct liners have tended to use GIT data for validation 
purposes. 8-11 Although these CAA models assume an anechoic (nonreflecting) termination, the suitability of this 
assumption within the context of the GIT has not been determined. For example, if the CAA and GIT results show' 
significant differences in the frequency range where the GIT termination is not anechoic, it would be difficult to 
determine whether such discrepancies w ere due to termination effects or other reasons. Educing the impedance with 
the measured and nonreflecting termination and demonstrating that the educed impedance is minimally affected could 
determine the suitability of using an anechoic boundary condition with GIT data. 

This paper provides a description of one of the first applications of the LEE model to data acquired in the modified 
GIT. GIT measurements obtained for ceramic tubular test liners arc used to provide a “sanity check” on the LEE 
model by showing that it reproduces the educed impedances of the more extensively tested CHE model in uniform 
flow. The current research also provides an opportunity for the authors to showcase some recently developed tools that 
significantly enhance the efficiency and accuracy of the impedance eduction methodology. Among these are a parallel, 
sparse, equation solver (that significantly reduces the wall clock time and computer memory requirements) and the 
Genetic algorithm (that adds increased robustness to the technique and renders it extremely attractive for massively 
parallel computing architectures). These new tools are used to extensively probe the acoustic solution in the vicinity 
of the duct wall and exphiin the anomaly in educed impedance at 0.5 kHz for high Mach number flows. To determine 
whether an anechoic termination is a suitable assumption for validating CAA codes using GIT data, the impedance 
educed using the measured termination is compared to that using an anechoic termination. 
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II. Statement of Problem 


F igure 1 is a schematic of the modified LaRC GIT and Cartesian coordinate system used in this study. The source 
and exit planes are located 203 mm upstream and downstream, respectively, of the leading and trailing edges of the 
test liner. The 406 imn long test liner (L2 — L\= 406 mm) is located on the upper wall of the GIT. The two sidewalls 
and bottom wall of the duct are rigid and the test liner is assumed uniform and locally reacting with an unknown 
dimensionless liner impedance, t. The cross section of the GIT is 51 mm in width and 51 mm in height. Although 
the duct exit phuie was designed to be anechoic, measurements show that reflections occur at sufficiently high mean 
flow Mach numbers. To account for these reflections, the exit plane impedance is measured and used as input to the 
impedance eduction method. It is assumed that the acoustic pressure field along the wall opposite the test liner is 
known from the wall mounted microphones (provided, for example, by the 95-microphone window 7 ). Finally there is 
a steady, nonunifonn, mean flow field that flows subsonically from left to right with mean flow Mach number, .Wo, as 
shown in Fig. 1. Ilie problem at hand is to determine the unknown impedance, £,, of the test liner in the presence of 
the flowing fluid. 


III. Duct Propagation Models 

I N this section a brief summary of two duct propagation models (the linearized Euler equations and convectcd 
Helmholtz equation models) that are solved to educe the liner impedance is presented. The physical geometry (uni- 
form cross section and test specimen, with rigid sidewalls) and restricted operational frequency range (/ = 0.5 — 3.0 
kHz) of the GIT, allow only two-dimensional acoustic disturbances to propagate. This situation is exploited by re- 
stricting the analyses to two space dimensions. Only a summary of the relevant differential equations and boundary 
conditions are presented. More detailed derivations may be found in any standard aeroacoustics text. 


A. Linearized Euler Equations (LEE) Model 

The governing differential equations for the 2-D LEE model are described by the 2-D linearized Euler equations 
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where the acoustic disturbance is assumed to take place isentropically in an ideal gas and a time dependence of the form 
c KII/ has been assumed. Implicit in the above derivation is that the mean flow field is fully developed (i.e.. Mo = Mq(x), 
po=constant, /b=constant, 7'o=constant). 

The proper forms of the boundary conditions have been well established. The wall impedance boundary condition 
for a locally reacting impenetrable surface of arbitrary shape and mean flow has been presented by Myers 12 
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where ^ = 00 along the rigid wall part of the GIT, and the unit normal vector to the duct wall is taken as positive when 
pointing into the fluid. The original form of the above boundary condition contains wall curvature terms. However, 
these curvature terms have been dropped because the walls in the GIT rue considered straight. At the source plane, the 
acoustic pressure and the transverse component of acoustic particle velocity are assumed known 


P = Ps(x), v = V s (x) 


( 5 ) 


In general, when sound propagating through the duct reaches the termination located at z = L, a radiated and a reflected 
wave field rue produced. This effect is best modeled for an assumed plane wave source by the exit impedance boundary 
condition 
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Finally, the Myers boundary condition 12 assumes the wall impedance is continuous along the liner. This assumption 
is violated in the GIT due to the abrupt impedance change at the liner leading and trailing edge (see Fig. 1). The LEE 
model therefore imposes continuity of the acoustic pressure and acoustic particle velocities at the leading and trailing 
edge of the liner to ensure that the acoustic pressure and velocity field remains a continuum. 


B. Convected Helmholtz Equation (CHE) Model 

In the CHE model the mean velocity field is taken as constant (Afo(.r)=constant) so that the acoustic pressure field in 
the GIT is governed by the convected Helmholtz equation 
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and the wall impedance boundary condition becomes 
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The sound source boundary condition is 

P = Ps 

and the exit impedance boundary condition becomes 

dp _ —ikp 
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( 9 ) 


( 10 ) 


Finally, the continuity of the acoustic pressure and its gradient at the leading and trailing edge of the liner is imposed to 
ensure that the acoustic pressure field remains a continuum. Equations (7)-(10) are simply a restatement of Eqs. (3)-(6) 
under the condition of uniform mean flow. 


IV. Solutions to the LEE and CHE Models 

A nalytical solutions for the sound field that satisfy the LEE and CHE models can be obtained only for certain 
simplified situations. Thus, one must resort to a numerical method that provides approximate, but acceptable 
solutions. The numerical method chosen to obtain the solution to both models is the conventional finite clement 
method. Both models employ a rectangular element ruid a uniform mesh. Specific details of the numerical method 
used in both models are described elsewhere. 3 5 However, to make this text as self-contained as possible, a brief 
summary of the numerical method for each model is presented. 

The convected Helmholtz equation model uses cubic Hermite polynomial basis functions and a weak formulation 
to satisfy the impedance boundary conditions. The LEE model uses linear polynomials for the basis functions and a 
strong formulation is used to satisfy the impedance boundary conditions. The sound source boundary condition for 
both models is implemented by constraining the nodal degrees of freedom at the source plane. 

The application of the finite element method to both the CHE and LEE model results in a system of simultaneous 
linear equations of the form 

[A]{d>} = {B} (11) 

Options were made available to solve equation (11) using a direct sparse solver for the CHE model and the LAPACK 
band solver for the LEE model (the sparse solver has not yet been implemented into the LEE model). Both solvers (the 
direct sparse and LAPACK band solver) run in parallel. The solution for the vector, {<!>}, is obtained on the computer 
platform Lomax. Lomax is an SGI ORIGIN' 3000 computer platform with 250 GB RAM and 512 CPUs. Each CPU 
runs at a speed of 600 MHz and uses the R 14000 processor chip. 

V. Impedance Eduction Technique 

T HE essential feature of the impedance eduction technique employed in this paper is the development of a general 
purpose numerical procedure to determine the impedance that reproduces the known acoustic pressure field on the 
wsill opposite the test liner. If the CHE and LEE model are indeed well-posed, then constructing a proof to show that 
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the wall impedance of the test specimen is indeed an extrema of the wall error pressure function is relatively straight 
forward. The procedure used here is to employ the wall error function 


nwall 



as an objective function and determine the unknown impedance such that this objective function is a minimum. The 
above objective function may be interpreted as the difference between the measured and finite element computed 
wall pressure opposite the test liner. Note that the proposed procedure is perfectly general and requires only that 
the CIIH and LEE formulation constitute a well-posed boundary v alue problem. Thus, nonuniform impedance and 
flows containing shear layers are possible within the context of the current impedance eduction theory. The reader 
is reminded, however, that the increased generality afforded by the current theory is partly mitigated by the need for 
a fast mid reliable equation solver (equation (1 1)) mid a robust method (i.e., optimizer) for minimizing the objective 
function. 


VI. Optimization Techniques 

I N this section three numerical strategies arc presented for determining the wall impedance that minimizes the ob- 
jective function. The impedance of the test liner is assumed constmit so that it consists of a two-parameter family, 
£ = 0-Mx- Thus, two parameters arc free to vary in the optimization: the dimensionless acoustic resistmice (0) mid 
dimensionless acoustic reactance (yj. 

A. Contour Deformation Method (CDM) 

The most reliable method for determining the minimum in the two-dimensional space defined by the uniform liner 
is the contour deformation method. In this method the dimensionless resistance is plotted along the abscissa mid the 
dimensionless reactance is plotted along the ordinate. Wall pressure error contours then represent the locus of points 
in the impedance plane yielding constant values of wall pressure error, F . These contours are then plotted using a 
graphing package. The impedance of the test liner is easily shown to be the location of the center of the contour 
with the smallest value of wall error. The CDM method is applicable only to a uniform impedance liner. It also has 
the disadvantage that it is not only computationally expensive, but also labor intensive because a computer-graphics 
program must be used to map the contours for each frequency mid mean flow Mach number. Therefore, two alternative 
numerical approaches for obtaining the minimum of the wall pressure error function are presented in this paper. These 
two alternative approaches me more automated and less labor intensive than the contour deformation method. They 
me also applicable to variable impedance liners. 

B. Stewart’s Adaptation of Davidon-Fleteher-Powell (SDFP) Method 

The first alternative approach to CDM is the Stewart’s adaptation of Davidon-Hetcher- Powell optimization method. 13 
The SDFP is representative of a general class of numerical algorithms known as descent methods. Its defining feature 
is that the impedance at the minimum point is obtained by iteration with the gradient of the objective being computed 
numerically by a finite difference approximation. It consists of a process for which the impedance is modified at 
each cycle of the iteration so that the objective function decreases in value. The iterative process is repeated until 
convergence to the minimum is achieved to within a specified tolerance. The SDFP combines the best features of 
steepest decent (good performance if mi initial value of the impedance is not close to the minimum point) and Newton’s 
method (rapid convergence if the initial value of impedance is close to the minimum point). The SDFP is generally 
considered to be a fast, reliable, optimization method and is expected to work well on the objective function used in 
this report. 

Although the convergence is usually rapid, the SDFP algorithm as applied in this paper has the disadvantages of 
requiring both an initial starting location for the impedmice mid the computation of the gradient of F. In general, the 
SDFP will locate only local minima. Thus, if the problem has multiple local minima, different initial values of the 
impedance may iteratively converge to different local minima. Also, when the minimum is located on a relatively flat 
surface in the function space, the educed impedance may vary because of differences in the convergence tolerance mid 
the direction from which the minimum is approach. Finally, the SDFP is a sequential algorithm that scales poorly on 
massively parallel computer architectures such as Iximax. 
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C. Genetic Algorithm (GA) 

The second alternative optimization algorithm to the CDM is the Genetic algorithm. 14 The GA is inspired by Darwin’s 
theory of evolution. Solution for the minimum point using the GA uses an evolutionary process. The algorithm begins 
with a randomly chosen set of liners called the initial population. The liners from the initial population me taken and 
used to form a new population of liners. Selection of the new population is motivated by the expectation that the new 
population will be better (i.e., have a smaller value of wall error) than the old population. Liners that are selected 
to form the new population (called the offsprings) are therefore selected according to their fitness (i.e., the value of 
F they produce)- the more suitable they are the more chances they have to reproduce. This reproductive process is 
repeated for a number of generations or until the condition, F is a minimum, is satisfied within some tolerance. 

In contrast to the SDFP, the GA does not require calculation of the gradient of the objective function, nor does it 
require an initial starting value for the impedance. Furthermore, the GA is a search method that searches the design 
space for a global minimal point. Thus, the principal drawback of the GA is that many more function evaluations are 
required when compared to that of SDFP. However, this is more than mitigated by the fact that the Genetic algorithm 
is highly scalable on massively parallel computing architectures and is normally run in parallel. A more complete 
description of both the SDFP and GA is beyond the scope of this paper. For a more comprehensive description of both 
automated search algorithms used in this paper, readers may consult Refs. 13 mid 14. 

VII. Test Liners and Mean Flow Conditions 

I N this paper impedances are educed for a CT65 and a CT57 ceramic tubular test liner. These two test liners are not 
typical of the cavity-backed, perforated face sheet construction presently used for noise abatement. They are chosen 
because their resistance is dominated by internal fluid viscosity mid is thus very insensitive to sound pressure level 
mid mean flow effects, mid their impedmice spectra vary over a range typical of liners used in aircraft engine nacelles. 
Figure 2 is a schematic of the CT65 mid CT57 ceramic liners. The CT65 liner consists of ’’sinusoid-shaped” parallel 
channels embedded in a ceramic matrix. These 77.5 mm-deep channels, with equivalent circular diameters of 0.76 
nun, mn perpendicular to the exposed surface and provide a surface porosity of 65 percent. The CT57 liner consists 
of a ceramic structure of parallel, cylindrical channels, embedded in a ceramic matrix. The channels are 85 mm deep 
(d = 85.0 mm) mid run perpendicular to the exposed surface. The diameter of the cylindrical channels, a , is chosen 
to provide a surface porosity of 57 percent. For both liners, each channel is rigidly terminated such that it is isolated 
from its neighbor to ensure a locally reacting structure. 

The impedance of each ceramic liner is educed at frequencies for which only the plane wave mode is cut on in the 
hard w all sections of the GIT. Results arc computed for source frequencies ranging from 0.5 to 3.0 kHz in increments 
of 0.5 kHz mid centerline Mach numbers, M c , of 0.000, 0.100, 0.200, 0.300, 0.475. A major complication in using 
GIT data is that in reality the flow contains a shear layer. Therefore, detailed mean flow profiles were acquired at three 
axial locations (z=203 imn, z=356 mm, mid '=559 imn) mid numerical integration w as used to compute mi average 
Mach number at each of the three selected axial planes. Finally, the average Mach number for each of the three axial 
planes w ere averaged together to obtained a “smeared” average Mach number, M ave , that is used as the uniform flow- 
Mach number for CHE mid LEE model. The average Mach number, dimensionless mean pressure, dimensionless 
mean temperature, mid dimensionless mean density measured in the GIT arc given in table 1 for each centerline Mach 
number, M c . The mean condition given in table 1 have been nondimensionalized with the constant mean conditions 
used in Ref. 15. Note that the measured values of mean pressure, memi temperature, and mean density varies somewhat 
with the centerline Mach number. This contrasts with the assumption used in Ref. 15 for which it was assumed that 
To = T re f — 295.0 K, Pq = P re f — 101325.0 Pa, and po = p re [ — L2 kg/m 3 . In the result section, a comparison is 
given of the impedmice educed using both the variable mean conditions (table 1) mid the constant mean conditions of 
Ref. 15. 


VIII. Optimization Parameters 

F OR each set of mean flow conditions the impedance is educed using both models (i.e., CHE and LEE model) mid 
all three optimization techniques (contour deformation (CDM), Stew art’s adaptation of Davidon-Fletcher- Pow ell 
(SDFP), mid the Genetic algorithm (GA)). Although a full description of all parameters used in each optimization 
method is beyond the scope of this paper a summary of salient features used in each method arc summarized below" 
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CDM Method 


1. A 201 x 101 evenly spaced impedance grid with an increment spacing of 0.05 (A0 = 0.05) on the dimensionless 
resistance axis and 0.1 (Ax = 0.1) on the dimensionless reactance axis 

2. Application of Tecplot (version 10) as the graphing package to determine the center of the contour containing 
the smallest value of F 


SDFP Method 

1 . No constraints on the dimensionless resistance and reactance 

2. A central difference approximation to the gradient of F 

3. All initial impedance of t = 2.0 — 1 .0/ 

4. A tolerance of 10 -8 as the stopping criteria 

GA Method 


1. Parallel function evaluations 

2. A population density of 100 liners for each generation 

3. An evolutionary process using tournament style selection and uniform crossover 

4. Creep and jump mutations 

5. An elitist strategy 

6. Constraints on the dimensionless resistance and reactance (0.0 < 0 < 10.0, — 10.0 < x < 10.0) 

7. Termination of the evolutionary process after 100 generations 

IX. Results and Discussion 

I N this section solutions obtained for the unknown impedance from the LEE model are compared to those of the 
more extensively tested CHE model for the range of frequencies and centerline flow Mach numbers for which 
measurements in the GIT could be obtained. Impedances are educed using each of the three optimization methods 
for an assumed plane wave source using the measured tennination impedance. One issue which one must come to 
grips with in this section is how to present the large number of results (two liners, three optimization techniques, four 
meshes, six frequencies, and six mean flow Mach numbers). Thus, 864 optimization runs were acquired anti only a 
selected number of these results can be presented. To better accommodate the massive number of optimization runs, 
results are presented in tabular instead of graphical formats. Graphical formats are presented only in a few selected 
cases for which a graph is appropriate. 

A. Parallel Solver Statistics 

Impedances educed for both the CHE and LEE models are obtained on a 65 x 11 baseline Unite element mesh 
(A z = 13 mm, mid Av = 5 mm). A mesh refinement study is also conducted on three additional meshes to ensure 
that the educed impedance has converged. The second mesh was obtained by doubling the density of the baseline 
mesh, the third by doubling the density of the second mesh, and fourth by doubling the density of the third mesh. 
Thus, the finest mesh has a density eight times finer than the baseline mesh. This subsection presents results from the 
part of the study concerning solver wall clock time (because Lomax contains 250 GB of RAM, memory requirements 
was not an issue). 

Table 2 shows the wall clock time in seconds, parallel speedup, and efficiency required to solve the finite element 
matrix equation (see Eq. (12)) for the sparse and LAPACK band solvers on the smallest mesh. Recall that the CHE 
model employs the parallel sparse solver, whereas, the LEE model employs the parallel LAPACK band solver. Al- 
though the matrix for (he sparse solver is 25 percent hirger (N = 1 66, 21 2) than that of the LAPACK solver (N = 124,659), 
the sparse solver obtains the solution using considerably less wall clock time. Note also that the parallel LAPACK 
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solver uses almost fifty percent more wall clock time on 64 CPUs than on 4. Similar statistics for the parallel speedup 
tmd parallel efficiency arc also given in the table. The reference timing for computing the parallel speedup is that 
required to solve the matrix equation on 4 CPUs (CPUs on Iximax me only allocated in multiples of 4). Note that 
parallel speedup is obtained for the LA PACK band solver on up to 16 CPUs, whereas, for the sparse solver the parallel 
speedup is obtained on at least 64 CPUs. Generally speaking, a parallel code is said to be efficient if it runs at 60 
percent efficiency or more. Applying this rule of thumb to statistics given in table 2, the parallel solvers mn efficiently 
on up to 8 CPUs. However, since both solvers show parallel speedup (i.e., reduction in wall clock) on as many as 
16 CPUs, it is reasonable to ran both solvers on as many as 16 CPUs to obtain speedup although they would not be 
running most efficiently. This strategy (running each solver on 16 CPUS) was employed in results obtained in this 
paper. 

In the foregoing paragraph, the statistics related to the parallel equation solvers were discussed. However, for a 
chosen mesh, the efficiency of the impedance eduction technique depends not only on the choice of solver but also 
upon the choice of a model. Results in this work show that the CHE model is considerably more efficient in educing 
the impedance on all four meshes than LEE. This result may reflect the fact that the CHE model has a higher order 
clement (cubic basis functions) than the LEE model (linear basis functions). For example, on the smallest mesh, the 
average wall clock time to educe the impedance (for each frequency Mid Mach number) using SDFP and the CHE 
model was 10.4 s. On the other hand, the average wall clock time required by SDFP using the LEE model was 17.4 s. 
Reduction in wall clock time on the smallest mesh was much more dramatic for the CHE when compared to that of 
the LEE model. The average wall clock time to educe the impedance (for a single frequency and Mach number) on the 
smallest mesh using SDFP with the CHE model was only 1,004 s, whereas, the LEE model required 10,228 s of wall 
clock time on this mesh. It was also observed that on the smallest mesh the CHE model consumed 1.17 GB of RAM 
compared to 1.73 GB for the LEE model. Such large reduction in wall clock time of the CHE model when compared 
to that of the LEE model cannot be explained by the differences in the efficiency of solvers alone, because each solver 
was run on 16 CPUs. 

B. Performance of Automated Optimizers 

Two different automated optimization algorithms (i.e., GA Mid SDFP) are used to educe the impedance of the test liner. 
The performance of the GA Mid SDFP on the baseline grid for the CHE model is now presented. Although results 
presented here Me for the CT57 liner at Mach 0.255, they arc representative of results obtained over the full rmige of 
Mach numbers Mid liners. Figure 3 shows the performance of SDFP for each frequency. Note that SDFP converges in 
seven to nine iterations with better performance at the higher frequencies. The w all error function at the anti -resonant 
frequency (2.0 kHz) is extremely fiat Mid SDFP is known to be deficient on flat surfaces. 13 Because SDFP use central 
differences to approximate the gradient of the objective function, approximately nine wall error function evaluations 
Me required per iterations. Thus, the time equivalent of 56 to 81 equation solves Me required before SDFP converges. 

The Genetic algorithm perfonnmice study is presented in Fig. 4. Recall that each population in the GA consists of 
100 offsprings. Therefore, the y axis value in Fig. 4 is the wall error corresponding to the offspring with the smallest 
wall error in the population. Note that the GA converges in ten generations or less Mid that progress in the minimization 
of the objective is extremely rapid over the first few generations. Further, the w all clock time required to compute a 
GA generation is approximately equal to that required to perform mi equation solve. It should be also be noted that 
just as w'ith SDFP, the performance of the GA is better at the higher frequencies (for example, the GA has converged 
in tw'o generations at 3.0 kHz). 

Generally speaking, in optimum design studies on pMallcl computing Mchitectures, the most importMit metric is 
the wall clock time required to obtain the optimum point. Thus, optimization algorithms that minimize the wall clock 
time Me normally preferred. Figures 3 Mid 4 provide insights as to which of the automated optimization algorithms 
returns the educed impedmice within the smallest wall clock time on Lomax. Note that the GA returns the optimum 
point in a wall clock time that is neMly mi order of magnitude less thmi SDFP. On the other hmid, the GA could not 
be run in pMallel on a sequential computing platform Mid would require 20% more wall clock time thmi the GA on a 
sequential computer platform. 

C. Infinitely Long CT65 Liner Results 

In the first exmnple, the LEE model is validated in a 5 1 -mmx5 1 -mm duct for mi infinitely long CT65 liner, for which 
the exact zero flow impedMice was known from normal incidence impedmice tube measurements. For simplicity, 
attention is focused only on downstremn propagation of the lowest order mode in the infinite duct. Source plane data, 
exit impedance, mid lower wall acoustic pressure data required to educe the impedance is synthesized from the exact 
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mode solution in an infinitely long duct. Figure 5 compares the dimensionless impedance educed using the CIIH model 
with the exact values using SDFP. Results arc plotted on a dmil axis system with resistance plotted along the yl axis 
mid reactance along the y2 axis. Because the CHE model educed the same impedance for each Mach number on each 
of the four meshes, only a single curve is shown. The numerically educed impedance from SDFP is in full agreement 
with the exact value. SDFP results using the LEE model were also computed on each mesh, mid results arc given in 
table 3. Good agreement between the exact mid educed impedances for each mesh is observed, with one exception. 
The baseline mesh (i.e., the coarsest mesh) has not converged to the correct impedance. This simply reflects the fact 
that the baseline mesh is too coarse for the impedance eduction process to be accurate. 

Before proceeding to the next example, it is worth noting that the Genetic algorithm mid the contour deformation 
method were applied to the infinitely long CT65 data. The educed impedances are not presented for the sake of 
brevity. However, the educed resistance and reactance from these two optimization algorithms were in agreement with 
that obtained by SDFP to within ±0.05. 

D. Finite Length CT57 Liner 

In this second example the impedance of the CT57 liner is educed using test data acquired at 140 dB reference sound 
pressure level in the GIT with the test liner installed. This example allows one to evaluate the models in the presence 
of real test data, with abrupt transitions at the leading mid trailing edge of the liner. Educed resistances for the CT57 
liner on the second mesh arc given in table 4, and the corresponding reactances arc given in table 5. The coarsest 
mesh results did not converge for the LEE model, but the second, third, and fourth mesh converged to nearly identical 
impedances. Therefore, only second mesh results me presented in the table (the CHE model converged for all meshes). 
Educed impedances me given for both the convectcd Helmholtz mid the linearized Euler models mid for each of the 
three optimization techniques. Results computed in Ref. 15 me also presented in the table. Recall that the Ref. 15 
results assumed the same set of memi conditions (7’o = T rc f , Pq = P re [ . and p 0 = p re f) for each memi flow Mach number. 
However, table 1 shows that these memi conditions change with the centerline Mach number. 

Several observations me noted from the CT57 educed impedances tabulated in tables 4 mid 5. Note that for each 
optimization method, the convected Helmholtz and Linearized Euler model educe essentially the same resistance mid 
reactance at each frequency mid flow Mach number with one exception. The two models educed different impedances 
at 0.5 kHz for high Mach numbers. This anomaly has been reported in previous papers and for different liners. It shows 
up here for the CT57 liner in the frequency and Mach number range (i.e., low frequency and high Mach number) for 
which shear layer effects are generally thought to be most important. Note also that the 0.5 kHz anomaly is not 
present at low Mach number and becomes more mid more dominant as the flow speed increases. Another important 
observation is that the results of Ref. 15 agree with results in this paper except where the anomaly occurs (i.e., 0.5 
kHz mid high Mach number) mid the anti-resonant peak (i.e., 2.0 kHz). 

Figures 6, 7, 8, mid 9 show contour maps of the reciprocal of the wall error function ( 1 / F) at 0.5 kHz for increasing 
flow speeds. The contour maps were produced from the convected Helmholtz model. Note that at zero Mach number 
(Fig. 6) the contour map has a unique, clearly defined eye that contains the maximum of 1 /F (i.e., the minimum of F). 
As Mach number increases (Fig. 7, 8, mid 9) the eye becomes more mid more diffused and the surface in the vicinity 
of the minimum is extremely flat with a dimneter of approximately 1 .2 poCo at the highest Mach number. Ilius, errors 
in educed dimensionless resistance and reactance of 1.2 W'ould not be unusual at this frequency. In all likelihood, the 
0.5 kHz contours at the high Mach numbers me non-physical mid probably generated by the absence of the shear layer 
that should be included in future computations. 

The effects of the duct termination on the educed resistance and reactance of the CT57 liner are given in table 6a 
mid 6b, respectively. Here the CHE model w'as used to evaluate the effects of the measured termination and for mi 
assumed anechoic termination using SDFP. Results in the table w ere computed on the second mesh. It is clem from 
these computations that both terminations give nearly identical educed impedances except at high Mach number near 
the shmp anti-resonant peak. Results similar to those of table 6 were found for the LEE model. 

Concluding Remarks 

This paper has presented a comparison of the impedances educed from a linearized Euler model to that of a well- 
tested convected Helmholtz model, under the assumption of a uniform memi flow field mid a plane wave source. 
Test data was acquired from ceramic tubular test samples installed in the NASA I .angles Grazing Incidence Tube. 
Impedances were educed using three search techniques (the contour deformation method, Stewart’s adaptation of the 
Davidon-Fletcher-Powell method, and the Genetic Algorithm) for the range of frequencies mid centerline Mean flow' 
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Mach numbers for which data consistent with the plane wave assumption could be acquired. The principal findings of 
this paper sue 

1 . All three optimization techniques consistently converged to the same impedance for all frequencies and average 
flow Mach numbers when used with the same model. 

2. Nearly identical impedances are obtained from the two models, with the exception that different impedances 
were educed at 0.5 kHz. at high mean flow speeds. 

3. The use of an anechoic termination with measured data in the NASA Langley Grazing Incidence Tube is valid 
except near a sharp anti -resonant peak. 

Consistency between the impedances educed using the two models also provides confidence that the linearized Euler 
model is ready for application to more realistic flows, such as those containing shear layers. Finally, because the 0.5 
kHz anomaly occurs at low frequency where the shear layers is known to be important, it is likely due to absence of 
the shear layer. 
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Tables and Figures 


Table 1. Measured Mean Flow Condition in the GIT 
(T ref = 295.0 K, Pref = 101325.0 Pa, and p re f = 1 -2 kg/m 3 ) 


M c 

Move 

wnusm 

wsusm 


0.000 

0.000 

1.003 

1.002 

1.001 

0.100 

0.079 

0.982 

1.006 

0.976 

0.200 

0.172 

1.021 

1.004 

1.017 

0.30 

0.255 

1.061 

1.005 

1.056 

0.400 

0.335 

1.130 

1.001 

1.129 

0.475 

0.400 

1.129 

1.004 

1.125 


Table 2. Sparse and LA PACK Solver Statistics for the Solution of [A] { <f> } = {/?} on the Finest Mesh 


Number 
of CPUs 

Sparse 
Wall Clock 
N= 166,212 

LAPACK 
Wall Clock 
N = 124,659 

Sparse 
Speedup 
N= 166,212 

LAPACK 
Speedup 
N = 124.659 

Sparse 
Efficiency 
N= 166,212 

LAPACK 
Efficiency 
N = 124,659 



74s 

1.00 

1.00 

100% 

100% 



60s 

1.29 

1.23 

65% 

62% 



65 s 

1.35 

1.14 

34% 

29% 


27 s 

77 s 

1.33 

0.96 

17% 

12% 


32 s 

108 s 

1.13 

0.68 

7% 

4% 
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Figure 2. Schematic of ceramic tubular test liners. 
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Table 3. SDFP educed impedance from LEE Model for CT65 liner 


Move 

/.kHz 

Exact 

Mesh 1 






0 

X 

e 

X 


X 


X 

0 

X 

0.000 

ana 

0.44 

-0.66 

0.45 

mm 

0.44 

-0.66 

0.44 

Esa 

0.44 

Era 

0.079 

M3 

0.44 

-0.66 

0.47 

Eg 

0.45 

-0.66 

0.45 

Eg 

0.44 

Eg 

0.172 

0.5 

0.44 

-0.66 

0.49 

-0.65 

0.46 

-0.66 

0.45 

-0.66 


-0.66 

0.255 

0.5 

0.44 

-0.66 

0.60 

-0.71 

0.48 

-0.69 

0.44 

-0.62 


-0.65 

0.335 

0.5 

0.44 

-0.66 

0.30 

-0.90 

0.56 

-0.72 

0.47 

-0.67 

0.45 

-0.66 

0.400 

0.5 

0.44 

-0.66 

0.52 

-0.51 

0.48 

-0.63 

0.46 

-0.65 

0.45 

-0.65 

0.000 

1.0 

0.40 

-0.61 

0.42 


0.40 


0.40 

ebi 

0.40 

EBI 

0.079 

1.0 

0.40 

-0.61 

0.45 

-0.59 

0.43 


0.41 

EH 

0.40 

ESI 


1.0 

0.40 

-0.61 

0.48 

-0.59 

0.43 

-0.61 

0.42 

-0.61 

0.41 

-0.61 


1.0 

0.40 

-0.61 

0.52 

-0.60 

0.45 

mm 

0.42 

-0.60 

0.41 

-0.61 

0.335 

1.0 

0.40 

-0.61 

0.63 

-0.52 

0.45 

ESI 

0.43 


0.40 

-0.61 

0.400 

1.0 

0.40 

-0.61 

0.47 

-0.47 

0.49 

-0.60 

0.43 


0.41 

-0.60 

0.000 

1.5 

0.63 

2.48 

0.67 

2.45 

0.64 

2.48 

0.63 


ESI 

2.48 


1.5 

0.63 

2.48 

0.64 

2.46 

Ml 

2.48 

0.63 


ESI 

2.48 


1.5 

0.63 

2.48 

0.59 

2.49 

Hi 

2.49 

0.61 

m 

EH 

2.48 


1.5 

0.63 

2.48 

0.61 

2.49 

0.59 

2.50 

EE1 

m 

0.61 

2.49 

0.335 

1.5 

0.63 

2.48 

0.57 

2.51 

0.57 

2.50 

M 

2.49 

0.61 

2.49 

0.400 

1.5 

0.63 

2.48 

0.53 

2.53 

0.55 

2.51 

EH 

2.49 

0.60 

2.48 

0.000 


0.73 

-1.47 

0.80 

-1.39 

0.74 

-1.46 

ESI 

-1.47 

0.73 

-1.47 



0.73 

-1.47 

0.89 

-1.39 

0.77 

-1.45 

m 

-1.46 

0.74 

-1.47 


2.0 

0.73 

-1.47 

0.97 

-1.36 

0.80 

-1.45 

Ksa 

-1.46 

0.75 

-1.46 


2.0 

0.73 

-1.47 

1.05 

-1.34 

0.83 

-1.44 

0.78 

-1.45 

0.76 

-1.46 

0.335 

2.0 

0.73 

-1.47 

1.17 

-1.26 

0.86 

-1.43 

0.79 

-1.45 

0.76 

-1.47 

0.400 

2.0 

0.73 

-1.47 

1.18 

-1.30 

1.04 

-1.59 

0.72 

-1.47 

0.75 

-1.48 

0.000 

2.5 

0.45 

0.17 

0.47 

0.13 

0.45 

EXE] 

0.45 


0.45 

0.17 

IEI 

2.5 

0.45 

0.17 

0.47 

0.15 

0.46 

Efl 

0.45 


0.45 

0.17 


2.5 

0.45 

0.17 

0.47 

0.17 

0.46 

0.18 

0.45 


0.45 

0.17 

m 

2.5 

0.45 

0.17 

0.47 

0.18 


0.18 

0.45 


0.45 

0.17 

0.335 

2.5 

0.45 

0.17 

0.47 

0.19 

ca 

-1.60 

0.45 

0.18 

0.45 

0.17 

0.400 

2.5 

0.45 

0.17 

0.49 

0.20 

m 

0.19 

0.45 

0.18 

0.45 

0.18 

0.000 

3.0 

0.49 

1.30 

0.45 

0.88 

0.54 

wm 

0.50 

wm 

0.49 

1.30 

0.079 

3.0 

0.49 

1.30 

0.52 

0.95 

0.53 

na 

0.49 

1 m 

0.49 

1.30 

0.172 

3.0 

0.49 

1.30 

0.67 

1.02 

0.35 

1.22 

0.48 

1.31 

0.48 

1.31 


3.0 

0.49 

1.30 

0.56 

1.09 

0.45 

1.32 

0.46 

1.31 

0.48 

1.31 

0.335 

3.0 

0.49 

1.30 

0.50 

1.03 

0.46 

1.32 

0.46 

1.32 

0.47 

1.31 

0.400 

3.0 

0.49 

1.30 

0.47 

1.04 

0.46 

1.33 

0.46 

1.32 

0.47 

1.31 
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Table 4. Educed Resistance for the CT57 Liner on the Second Mesh 


Move 

/.kHz 

Convccted Helmholtz 

Linearized Euler 



GA 

SDEP 

CDM 

Ref. 15 

GA 

SDFP 

CDM 

0.000 


0.55 

0.53 

0.53 

0.58 

0.55 

0.53 

0.53 

0.000 


0.47 

0.49 

0.48 

0.49 

0.47 

0.49 

0.49 

0.000 

1.5 

1.18 

1.16 

1.15 

1.12 

1.18 

1.16 

1.15 

0.000 


5.12 

5.12 

5.10 

4.75 

5.04 

5.08 

5.09 

0.000 


1.50 

1.47 

1.48 

1.49 

1.50 

1.47 

1.45 

0.000 

3.0 

0.71 

0.75 

0.74 

0.74 

0.71 

0.75 

0.75 

0.079 

0.5 

0.71 

0.71 

0.71 

0.77 

0.63 

0.63 

0.66 

0.079 

1.0 

0.39 

0.40 

0.41 

0.39 

0.47 

0.44 

0.42 

0.079 

1.5 


1.08 


1.00 

1.10 


1.07 

0.079 

2.0 


4.10 

■H 

3.69 

4.09 


4.11 

0.079 

2.5 


1.38 

1.37 

1.45 

1.42 

1.38 

1.39 

0.079 

3.0 


0.70 

0.71 

0.68 

0.71 

0.72 

0.71 

ESI 

0.5 

ESI 

0.79 

mm 

0.85 

0.63 

0.63 


m 

1.0 

m 

0.33 

B 

0.33 

0.39 

0.41 


0.172 

1.5 

1.26 

1.26 


1.22 

1.26 


1.25 

0.172 

2.0 

4.96 

4.94 

sa 

4.75 

4.96 


4.94 

0.172 

2.5 

1.10 

1.11 

1.10 

1.15 

1.10 

1.12 

1.11 

0.172 

3.0 

0.79 

0.81 

0.80 

0.79 

0.87 

0.83 

0.85 


0.5 

0.63 

0.62 

0.63 

0.64 

0.47 

0.45 

0.47 

0.255 

1.0 

0.24 

0.27 

0.26 

0.27 

0.39 

0.39 

0.40 

0.255 

1.5 

1.26 

1.24 

1.24 

1.18 

1.26 

1.23 

1.22 

0.255 

2.0 

5.35 

5.33 

5.32 

5.31 

5.35 

5.35 

5.33 

0.255 

2.5 

0.94 

0.97 

0.99 

1.02 

0.94 

0.98 

1.00 

0.255 

3.0 

0.79 

0.76 

0.75 

0.67 

0.79 

0.78 

0.80 


0.5 

0.71 

0.69 

0.70 

0.70 

0.32 

!KQ 

0.31 


1.0 

0.24 

0.20 

0.20 

0.20 

0.39 



0.335 

1.5 

1.18 

1.19 

1.20 

1.19 

1.18 

tg 


0.335 

2.0 

4.57 

4.60 

4.59 

4.81 

4.17 

B 

4.19 

0.335 

2.5 

0.94 

0.91 

0.90 

0.92 

0.87 

B 

0.93 

0.335 

3.0 

0.71 

0.73 

0.72 

0.72 

0.79 

S£3 

0.78 

0.400 

0.5 

0.87 

0.86 

0.88 

0.88 

0.47 

0.01 

0.47 

0.400 

1.0 

0.08 

0.07 

0.05 

0.07 

0.32 

0.30 

0.29 

0.400 

1.5 

1.26 

1.26 

1.27 

1.23 

1.34 

1.32 

1.30 

0.400 

2.0 

2.44 

2.47 

2.45 

2.77 

2.44 

2.41 

2.40 

0.400 

2.5 

0.79 

0.80 

0.79 

0.84 

0.79 

0.78 

0.79 

0.400 

3.0 

0.79 

0.75 

0.76 

0.66 

0.79 

0.78 

0.79 
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Table 5. Educed Reactance for the CT57 Liner on the Second Mesh 


Move 

/.kHz. 

Convected Helmholtz 

Linearized Lulcr 



GA 

SDFP 

CDM 

Ref. 15 

GA 

SDFP 

CDM 

0.000 



-1.70 

-1.70 

-1.74 

-1.69 

-1.70 

-1.69 

0.000 



0.13 

0.13 

0.13 

0.12 

0.13 

0.13 

0.000 

1.5 


1.28 

1.27 

1.29 

1.27 

1.28 

1.27 

0.000 


0.90 

0.88 

0.88 

0.81 

0.96 

0.85 

0.84 

0.000 


-1.61 

-1.59 

-1.58 

-1.59 

-1.61 

-1.58 

-1.59 

0.000 

3.0 

-0.20 

-0.18 

-0.19 

-0.19 

-0.20 

-0.19 

-0.18 

0.079 


-1.58 

-1.55 

-1.55 

-1.58 

-1.61 

-1.59 

-1.58 

0.079 


0.04 

0.08 

0.09 

0.08 

0.04 

0.07 

0.08 

0.079 

1.5 

1.29 

1.27 

1.28 

1.26 

1.22 

1.24 

1.22 

0.079 


0.98 

1.00 

0.98 

1.17 

0.90 

0.92 

0.91 

0.079 


-1.45 

-1.48 

-1.47 

-1.45 

-1.45 

-1.47 

-1.46 

0.079 

3.0 

-0.27 

-0.26 

-0.26 

-0.28 

-0.27 

-0.26 

-0.27 

EBa 


-1.50 

-1.49 

-1.49 

-1.50 

-1.61 



Eta 


0.20 

0.17 

0.16 

0.17 

0.12 



0.172 

1.5 

1.29 

1.31 

1.30 

1.31 

1.22 


1.24 

0.172 


-0.51 

-0.54 

-0.55 

-0.21 

-0.67 

-0.67 

-0.68 

0.172 


-1.45 

-1.48 

-1.49 

-1.47 

-1.45 

-1.49 

-1.48 

0.172 

3.0 

-0.12 

-0.14 

-0.15 

-0.15 

-0.20 

-0.16 

-0.17 

mm 

0.5 

-1.06 

-1.07 

-1.05 



mm 

■SI 

0.255 

1.0 

0.12 

0.11 

0.10 





0.255 

1.5 

1.37 

1.38 

1.39 

1.37 

1.29 

1.32 

H9 

0.255 

2.0 

-1.37 

-1.34 

-1.34 

-1.05 

-1.45 

-1.45 

-1.47 

0.255 

2.5 

-1.45 

-1.44 

-1.45 

-1.44 

-1.45 

-1.48 

-1.50 

0.255 

3.0 

-0.20 

-0.20 

-0.20 

-0.12 

-0.20 

-0.21 

-0.20 


0.5 




-0.57 

-0.98 

-1.31 

-1.30 


1.0 




0.15 

0.04 

0.07 

0.08 

0.335 

1.5 

1.29 

1.32 

1.32 

1.31 

1.29 

1.30 

1.29 

0.335 

2.0 

-1.69 

-1.67 

-1.67 

-1.62 

-1.69 

-1.70 

-1.69 

0.335 

2.5 

-1.37 

-1.39 

-1.40 

-1.39 

-1.45 

-1.46 

-1.46 

0.335 

3.0 

-0.20 

-0.23 

-0.24 

-0.23 

-0.27 

-0.24 

0.24 

0.400 

0.5 

-0.35 

-0.34 

-0.35 

-0.35 

-0.75 

-1.01 

-1.00 

0.400 

1.0 

0.12 

0.14 

0.12 

0.14 

0.04 

0.05 

0.07 

0.400 

1.5 

1.14 

1.16 

1.13 

1.19 

1.22 

1.19 

1.20 

0.400 

2.0 

-2.63 

-2.63 

-2.65 

-2.63 

-2.47 

-2.43 

-2.40 

0.400 

2.5 

-1.29 

-1.25 

-1.23 

-1.25 

-1.29 

-1.27 

-1.26 

0.400 

3.0 

-0.27 

-0.24 

-0.23 

-0.23 

-0.27 

-0.28 

-0.27 
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Tabic 6a. Educed Resistance for the CT57 liner using the Measured mid Anechoic Termination 


f, kHz 

Mo = 

0.000 

M 0 = 

0.079 

Mo = 

0.172 

Mo = 

0.255 

M> = 

0.355 

Mo = 

0.400 


Mcas 

PcK’o 

Mcas 

PoO> 

Meas 

Pof'o 

Meas 

PoO> 

Meas 

P()<"0 

Meas 


0.5 

0.53 

0.51 

0.63 

0.62 

0.63 

0.65 

0.45 

0.43 

0.35 

Eg 

0.43 

0.43 

1.0 

0.49 


0.44 

0.44 

0.41 

0.41 

0.39 

0.39 

0.36 

m 

0.30 

0.30 

1.5 

1.16 

1.15 


1.08 

1.27 

1.27 

1.23 

1.23 

1.17 

m 

1.32 

1.35 

2.0 

5.08 

4.97 


4.08 

4.95 

4.98 

5.35 

5.30 


H 

2.40 

2.45 

2.5 

1.47 

1.46 

1.38 

1.38 

1.12 

1.12 

0.98 

0.98 

B 

0.90 

0.78 

0.78 

3.0 

0.75 

0.75 

0.72 

0.72 

0.83 

0.83 

0.78 

0.78 

B a 

0.77 

0.78 

0.78 


Table 6b. Educed Reactance for the CT57 Liner using the Measured and Anechoic Termination 


fjcllz 

M 0 = 

0.000 

M 0 = 

0.079 

M 0 = 

0.172 

Mo = 

0.255 

M 0 = 

0.355 

M 0 = 

0.400 


Meas 

Po<~o 

Meas 

Po<"o 

Meas 

Po<"o 

Meas 

PcKo 

Meas 

Pof'o 

Meas 


0.5 

-1.70 

-1.69 

-1.59 

-1.59 



-1.35 

-1.45 

-0.96 

-0.95 

-0.76 


1.0 


0.13 

0.07 

0.07 



0.07 

0.07 

0.07 

0.07 

0.05 

0.05 

1.5 

1.28 


1.24 

1.26 



1.32 

1.32 

1.30 

1.29 

1.19 

1.22 

2.0 

0.85 


0.92 

0.97 

-0.67 

-0.75 

-1.45 

-1.79 

-1.70 

-1.90 

-2.43 

-2.49 

2.5 

-1.58 

-1.58 

-1.47 

-1.47 

-1.49 

-1.49 

-1.48 

-1.48 

-1.46 

-1.46 

-1.27 

-1.27 

3.0 

-0.19 

-0.19 

-0.26 

-0.26 

-0.16 

-0.16 

-0.21 

-0.21 

-0.24 

-0.24 

-0.28 

-0.28 
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Wall error function 



Figure 3. Performance of SDFP algorithm for CT57 test liner oil baseline grid (,\/o = 0.255). 
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Wall error function 



Figure 4. Performance of GA algorithm for CT57 test liner on baseline grid ( Wo = 0.255). 
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Figure 5. SDFP educed impedance from (he Conceded Helmholtz model for (he CT65 liner. 
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Dimensionless reactance 



Dimensionless resistance 
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Figure 6. Contour Map of the wall error function at 0.5 kHz for the CT57 liner (A/o = 0.0). 


21 of 24 



Dimensionless Reactance 
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Figure 7. Contour Map of the wall error function at 0.5 kllz for the CT57 liner (Mo = 0.079). 
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Dimensionless reactance 



F 

5.6763 

5.2998 

4.9233 

4.5467 

4.1702 

3.7937 

3.4172 

3.0406 

2.6641 

2.2876 

1.9110 

1.5345 

1.1580 

0.7814 

0.4049 


Dimensionless resistance 


Figure 8. Contour Map of the wall error function at 0.5 kHz for the CT57 liner (Mq = 0.255). 
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Dimensionless reactance 
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Figure 9. Contour Map of the wall error function at 0.5 kHz for the CT57 liner (A/q = 0.400). 
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